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ABSTRACT 

As a step toward a complete theoretical integration of 3D compressible hydrody- 
namic simulations into stellar evolution, convection at the surface and sub-surface layers 
of the Sun is re-examined, from a restricted point of view, in the language of mixing- 
length theory (MLT) . Requiring that MLT use a hydrodynamically realistic dissipation 
length gives a new constraint on solar models. While the stellar structure which results 
is similar to that obtained bv YREC (iGuenther. Demaraue. Kim &: Pinsonneault I1992I : 
Bahcall Sz Pinsonneault 20041 ) and Garching models (jSchlattl. Weiss fc Ludwiglll997l l. 
the theoretical picture differs. A new quantitative connection is made between macro- 
turbulence, micro-turbulence, and the convective velocity scale at the photosphere, 
which has finite values. The "geometric parameter" in MLT is found to correspond 
more reasonably with th e size of the strong downward plumes which drive convection 
( Stein Sz Nordlund 1998 1). and thus has a physical interpretation even in MLT. Use of 
3D simulations of both adiabatic convection and stellar atmospheres will allow the de- 
termination of the dissipation length and the geometric parameter (i.e., the entropy 
jump), with no astronomical calibration. 

A physically realistic treatment of convection in stellar evolution will require addi- 
tional modifications beyond MLT, including effec ts of kinetic energy flux, entrainment 
(the most dramatic difference from MLT found by Meakin Sz Arnett ( 2007bl ) ), rotation, 
and magnetic fields ( Balbus fc Hawleylll998l : iBalbus 2008 ). 



Subject headings: stars: evolution - stars: hydrodynamics - convection - Sun: photosphere 
white dwarfs - atmosphcres:3D - Binaries: eclipsing 



1. Introduction 

Recent simulations of three-dimensional 
compres sible convection and their theoretical 
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This in ixing length £ is parameterized as 



mixing-length theory (MLT), as currently 



^Steward Observatory, University of Arizona, 933 
N. Cherry Avenue, Tucson AZ 85721 

^School of Earth and Space Exploration, Arizona 
State University, Tempe, AZ 



is the local pres- 
sure scale height, and uml is adjusted to 
reproduce the radius of the present-day Sun. 
However, instead of being an adjustable pa- 
rameter, the mixing length is found to cor- 
respond to the dissipation length of the 
turbulence ( Kolmoeorov Il94ll . ll962l ). and de- 
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termined by the size o f the largest eddies 

dMeakin k Arnettll2007bl : lArnett. Meakin. Younel 
200S ) . Prom our own simulations (jMeakin &: Arnett 
2003 ) and those of others we find a robust ten- 
dency for the dissipation length to be 



numerical simulations and by laboratory ex- 
periment (e.g., fluctuations, non-locality, and 
entrainment), the framework of the standard 



id w min(£cz,4i?p). 



(1) 



where Icz is the depth of the convection zone. 
For shallow convection zones, the dissipation 
length is limited by the depth of the convec- 
tive region, and seems to approach a limiting 
value of id « AHp for deep convection zones. 

If aML is fixed, other parameters in MLT, 
which are generally left fixed by histori- 
cal convention, may be adjusted to com- 
pens a te fe.g.. (ITassoul. Fontaine. &: Winget 



19901 : ISalaris k Cassisil 120081 ^'). The most sig- 
nificant of these parameters is the geometric 
factoi0, which adjusts the rate at which radia- 
tion limits the degree of entropy excess in the 
super-adiabatic region (SAR). For simplicity 
we use qml to denote the geometric factor in 
units of the value used in conventional MLT 
(see Appendix for details). If we identify the 
geometric parameter as a measure of the size 
of the SAR, we remove the last free parame- 
ter in MLT. Although MLT is an incomplete 
theory, it does serve as a useful "language" 
to explain some of the changes implied by 3D 
simulations. 



The mixing length theo ry itself (jVitense 
1953 : Bohm-Vitense 19581 ). if used consis- 
tently, does capture many (but not all) as- 
pects of turbulent convection. However, a real 
replacement for MLT will provide a global so- 
lution and relax the local connection between 
the superadiabatic gradient and the enthalpy 
flux, so that regions of the convective zone can 
be subadiabatic, as observed in simulations. 

In order to establish a "baseline" from 
which to compare new effects demanded by 



^This is essentially t he c factor 
iTassoul. Fontaine, fc Wingetl (|l990D . 



solar model (jBahcall. Serenelli. k Pinsonneault 
20041 ) is examined with respect to modifica- 
tion of some aspects of convection. In this 
paper we show that the addition of dynami- 
cally realistic values of mixing length and ge- 
ometric factor give some interesting insights 
into the nature of the average stratification 
of the Sun just below the photosphere. In 
Section 2 we construct a series of solar mod- 
els with oml — Oml pairs, to delineate their 
properties. The notion of a "standard so- 
lar model" derives from the work of John 
Bahcall and col l abora tors, and is summa- 
rized in iBahcalll (|l989l l. It represents what 
is probably the most carefully tested aspect 
of the theory of stellar evolution. In Section 3 
we compare our models to standard results 
using the Yale Rotational Evolution Cod e 
(iGuenther. Demargue. Kim Pinsonneault 

19921: |BahcalLfcPinsonneault''2004l ), and the 

Garch ing code (see|Schlattl (1996); Schlattl. Weiss k Ludwie 
(1993)). In Section 4 we compare the outer 
layers of our models t o the 3D atmosphere s 
of Nordlund and Stein (jAsolund. et al. II2005I '). 
semi-empirical ruodels of the solar atmosphere 



(jPontenla et al. I l2006l l. and re-examine the 



question of convective velocities at the photo- 
sphere. In Section 5 we summarize the impli- 
cations of this work. 

2. Solar Models 

2.1. Hydrodynamically-consistent MLT 
Parameters 



Tassoul. Fontaine, k Wingetl (jl990l l have 



defined the parameters in MLT in a concise 
way: they define three parameters a, b, and 
c in terms of an a djustable mixing length i. 
In the notation of lArnett. Meakin. k Young 
(|2009l ). AV = V - Ve, so we have 



of 



a = vlHp/fgpT^V, 



(2) 
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Table 1 
Mixing Length Parameters^ 



MLT choice 


a 


b 


c 






0.125 


0.5 


24 


free parameter 


ML2= 


1 


2 


16 


free parameter 


AMY<i 


ss 0.1 


0.256 


24/.9ML° 


min(4iJp,£c'z) 



^As defined in lTassoul. Fontaine, fc Winged (|l990l ). 
^B6hm-Viteiis3 (|l958h : this is MLl. 

TassouL Fontaine, fc Wingetl (|l990| ): ISalaris fc Cassisil (|2008fl 



lArnett. Meakin. fc Youii3 (|2009t ). and this paper. 



'^guL = i^/V^Tb)"^, where rf, is the radius of a blob just contained 
inside the SAR. 



F,Hp/pvcCpT£AV, 



(3) 



and 



C = Cpp'^£VcK{Ve - Va)/(TT^AV (4) 



Table [T] gives standard values for the mix- 
ing length parameters in th e formu lation of 
Tassoul. Fontaine, Wingetl (|l990l l. The 
first t wo entries are t he st andard " flavor" 
due to lBohm-Viteiis3 (Il958l). and th e ML2 
flavor of Tassoul. Fontaine, Winget ( IQDOl ). 
In both cases the values of a, b, and c are 
fixed and the mixing length £ adjusted to 
reproduce the solar radius of the present 
day Sun. The third line presents the val- 
ues of these para meters as estimat e d froin 
3D simulations dMeakin fc ArnettI l2007bl : 
Arnett. Meakin. fc Young] I2OO9I ). Two strik- 
ing differences are apparent: (1) the mix- 
ing length is not an arbitrary constant. For 
deep convection zones (like the Sun which 
has icz = 20 Hp) the mixing length (dis- 
sipation length) approaches AH p. (2) the 
"c" parameter is intimately related to the 
geometric factor, and we assume to the 
thickness of the superadiabatic layer. Both 
the "a" and "b" parameter are fixed at 



the values for adiabatic turbulent convec- 
tion (jArnett. Meakin. fc Youn^l2009l ) . leaving 
c = 24/ gML as the remaining free parameter. 
Note that the flavors BV58 (MLl) and ML2 
both differ from those suggested by the simu- 
lations. 

Even though there are additional effects 
shown in 3D turbulent simulations which are 
not contained in MLT, it is useful to examine 
those changes which can be captured with a 
standard stellar evolution code. If the value 
of oml = ^/Hp is fixed, which of these pa- 
rameters is to be varied in MLT to get an 
acceptable solar model? The only parame- 
ter sufficient to the task is the "geometric 
parameter" (i.e., c or qml)- A simple way 
to examine the effects of the geometric pa- 
rameter is to vary its value relative to the 
value used in conventional MLT; we denote 
this scaled value by qml (see the Appendix 
below). We relate this factor to the size of the 
SAR by QML = (4/\/3rfe)2, .^j^^^.^ ^^le "blob 
diameter" 2rb is the thickness of the SAR. 
MLT results if we set ^(a/l = 1; this identifies 
th e SAR with the supera d iabati c "element" 
of Kippenhahn fc Weigert ( 199d ). p. 50. In 



3 



MLT, the "blob" is assumed to have a di- 
mension fixed by the mixing length. This is 
inconsis tent with 30 atm osphere simulations 
( Stein ^Nordlundl Il998l)') and solar model s 



Guenther. Demargue. Kim &: Pinsonneault 



1992 ): ISchlattl. Weiss Ludwid (|l997l l and 
below), which show that the superadiabatic 
region is narrow, less than a pressure scale 
height thick. MLT has two characteristic 
lengths, one of which is ignored by forcing 
the geometric parameter length scale to be 
the same as the mixing length (turbulent dis- 
sipation length). We will allow the "blob" 
size to differ from the mixing length in order 
to vary guL- 

For theoretical clarity we will apply ra- 
diative diffusion theory consistently up to 
the photosphere. While radiation transfer 
theory is more accurate than radiative dif- 
fusion, it is more cumbersome, and itself 
i s affected by the conv ection model used 
( VandenBerg. et al. 20081 ). After we under- 
stand the convection problem better, this ap- 
proach can be extended by a more sophisti- 
cated multi-dimensional treatment of radia- 
tive transfer in the outer regions. Compar- 
ison to 3D hydrodynami c atmospheres (e.g., 
Stein fc Nordlundl (Il998l 'l: iNordlund fc Stein 
ted)) can test the validity of this approach. 



2.2. Standard Input Physics 

These computations were done using the 
TYCHO stellar evolution code (revision 12; 

version control by SVN). Op acities were from 

Iglesias fc Rogers! (119961) andlAlexander Sz Ferguson 



(Il994l) with iGrevesse fc Sauvall (Il998ll abun 



(|l99d ). with the modifications via the di- 
mensionless geometric factor qml as given 
in the Appendix; qml = 1 gives conven- 
tional MLT. We stress that convection is 
treated in exactly the same way, with the 
same par ameters, in the interior a nd in the 
envelope ( VandenBerg. et al. 20081 ). Diffu- 
sion was treated with the Thoul subrou- 



Thoul. Bahcall. k Loeb 


1994)) 


ievitation ( 


Michaud et al. 




2004| 



radiative 

was ignored. The nuclear reactions were 
solved in a 177 isotope network using Reaclib 
(iRauscher fc Thielemanl]l200d ): weak screen- 
ing rates were incorporated as in John Bah- 
call's exportenergy.f program. The changes 
in metalicity due to nuclear reactions and to 
diffusion were taken into account by interpo- 
lation in both the opacity and the equation of 
state tables. The same equation of state and 
opacity tables are used in the interior and the 
atmosphere. 

2.3. Modified MLT Models 

We examine the solar models resulting 
from several choices of the mixing length, each 
constructed by varying the geometric factor 
qm l until the correct radius of the present day 
Sun was obtained. The other MLT parame- 
ters a and b are the flavor MLl in Table [TJ 

Five such models were constructed, with 
values of uml ranging from 1.6 to 5.2, as 
summarized in Table [5J The model A has 
OiML = 1-643, Qml = 1-0, and is typical 
of current solar models which use the Ed- 
dington gray atmosphere as the outer bound- 
ary condition and conve n tiona, l MLT (e.g., 



danc es. The OPAL-EOS (jRodgers. Swenson Iglesijfilchlattl. Weiss fc Ludwigl (|l997l )). This pro 
19961 ) equation of state was used over the 
range of conditio n s rel evant here. The 
Timmes fc Swestv (I2OOOI ) equation of state 
is automatically used for higher densities 
and temperatures, with a smooth interpola- 
tion across the joinin g region. The formula- 
tion of MLT is from iKippenhahn &: Weigert 



vides a baseline for comparison. Model D has 
OiML = 4.0, which is most consistent with hy- 
drodynamic simulations. 

Arnett. Meakin. fc Youn3 (j2009l ) found 



that aML was not constant, but depended 
upon the flow properties, and the equation 
of state. For solar models the surface con- 
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Table 2 
Solar models with MLT 



Model 


a.ML 


9ml 


R/Rq 


L/Lq 


rcz/Re 


HCsurf 




A 


1.650 


1.0 


1.001 


1.000 


0.7169 


0.2379 


2.25 


B 


2.323 


42.0 


1.001 


1.000 


0.7172 


0.2378 


2.80 


C 


3.286 


270.0 


1.001 


0.9997 


0.7173 


0.2377 


3.05 


D 


4.000 


595.0 


1.000 


0.9997 


0.7168 


0.2373 


3.20 


E 


5.190 


1,540.0 


1.001 


0.9998 


0.7172 


0.2377 


3.40 


Sun 






1.000 


1.000 


0.713 ±0.001 


0.24 


3.20^^ 



'■Inferred from the model data in lAsplund. et al. I (|2005r ). 



vection zone is deep, and changes little, so 
taking a constant oml is an adequate ap- 
proximation for this particular example. In 
MLT, the velocity obtained by a convective 
eddy is computed from the work done by the 
buoyancy force over a mixing length: 



aljL9Hp(3TiAV)/8, 



(5) 



where Pt is the compressibility. Hp the pres- 
sure scale height, and AV = V — Vad is 
the usual "super-adiabatic excess." For a 
given convective luminosity, larger uml im- 
plies larger velocities. 

Shallow convection zones, having shorter 
distances for buoyant acceleration to work, 
will have smal ler values of a m and smaller 
velocity scales ( Arnett. Meakin. fc Younglbood ) 
As the depth of the convection zone increases, 
the size of the largest eddies also rises, imply- 
ing larger uml- Such an increase will not con- 
tinue indefinitely; more vigorous convection 
develops more violent dissipation. The value 
of uml seems to "saturate" for very deep 
conve c tion zones (lArnett. Meakin. &: Young 
20091 : iMeakin k ArnettI booof h The solar 
convection zone is 20 pressure scale heights 
deep, and has yet to be simulated for its 
full depth with resolutio n a s high as used in 
Meakin fc ArnettI toOJlh or lStein fc Nordlund 
( 19981 ). Here we will examine the case in 



which such saturation occurs at oml ~ 4. 
This may be appropriate for the simulations 
of Nordlund and Stein (R. S tein, private com- 
munic ation) and those of IMeakin &: Arnett 
( 20091 ), and is consistent w i th the insensitivity 
of the Stein Sz Nordlund ( 19981 ) simulations 
to the exact position of the lower boundary, 
which was deeper than this . Further analysis 
of thi s issue is in progress IMeakin &: Arnett 
(|2009l l: 3D simulations for convective zones 
of depth 0.5 to 5 pressure scales heights seem 
consistent with this interpretation. The dis- 
tribution of values for om l in Table [2] covers 
this range. 

Softer equations of state, such as in par- 
tial ionization zones or electron-positron pair 

zones, give less vigorous veloci ties, but do not 

change the qualitative pic ture (I Arnett. Meakin. &: Young 
20091) . The simulations of lPorter &: Woodward 



(2000), : or an ideal gas equation of state, also 
seem to suggest that saturation may be be- 
ginning around uml ~ 3, which is consistent. 

Figure [T] shows the evolutionary tracks in 
the HR diagram, for each of the models (which 
differ only by the uml and qml parame- 
ters). For each aML-, a value of qml is cho- 
sen which gives a reasonable radius for the 
present-day Sun. After passing the stellar 
birthline (deuterium burning, L 10 Lq, see 
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Stabler Sz Fallal (|200J)), the tracks are very 
similar for all five models. However the in- 
creasing values of uml imply increasing tur- 
bulent velocities (Eq. [5]) . Table [2] gives values 
of the radius (R/Rq), luminosity (L/Lq), the 
radius of the lower boundary of the convec- 
tive zone (rcz/Re), the surface (convective 
zone) abundance of helium by mass fraction 
(Hesurf)) and the maximum mean turbulent 
velocity in the convection zone {vm) in km/s. 
The increase in Vm with ckml is clear. 

The models were adjusted to radius and 
luminosity of the present-day Sun to about 
one part in 10^ or better, which is sufficient 
to show accurately the differential effects to 
be discussed here. 

It is well known that, once a calibration 
of MLT parameters is done to fit the so- 
lar radius, paths in the HR diagram are 
little affected by which parameters were 
used (Federsen. Vandenberg. Sz Irwin 1990l : 



Salaris Sz Cassisi 



200 



However, the vari- 
ati on of the velocity scale, although not iced 
by iFedersen. Vandenberg. &: IrwinI (|l990i ) for 
example, has not been stressed. In Table [2] it 
is striking that only the velocity scale varies 
significantly with the variation of oml and 
9ml pairs constrained to fit the radius and 
luminosity of the present day Sun. This ve- 
locity scale is crucial for rates of entrainment, 
wave generation, and mass loss, and so may 
ultimately cause a change in the evolution- 
ary behavior when such effects are correctly 
included. It will be argued below that this 
variation in the velocity scale has direct ob- 
servational consequences (via line profiles and 
micro- and macro-turbulent velocities). 

The values of the lower radius of the solar 
convection zone rcz and the surface helium 
abundance Hcsurf are slightly different from 
the values of the standard solar model. Fart 
of the difference may be due to small errors in 
our stellar evolution code, which is not yet in 
its fully verified state. However the standard 



solar model uses MLT and therefore ignores 
several significant aspects of convection: tur- 
bulent heating, flux of kinetic energy, and en- 
trainment. These effects may move our mod- 
els toward the inferred values from helioseis- 
mology. In any case the differential effects we 
discuss here are much larger, and unlikely to 
be affected by small modifications in the refer- 
ence model. Notice that the predicted values 
of rcz and Hcsurf vary in only the fourth sig- 
nificant figure for models A through E while 
the velocity scale increases by more than 40 
percent. 

3. YREC and Garching Solar Models 

Our solar models are in good agreement 
with YREC and Garching models, but not 
yet as close to either as they are to each 
otheiH In this paper we concentrate on the 
differences caused by changing MLT param- 
eters and outer boundary conditions, rather 
than finding the absolute best standard so- 
lar model. Finding the absolute best solar 
model requires going beyond present formula- 
tions used in YREC and Garching codes (e.g., 
Michaud et"all (l2004l):lMeakin &: ArnettI (l2007bll: 
Arnett. Meakin. Sz Yound (|2009l l): we plan to 
address this in detail in future publications. 

3.1. Empirical Outer Boundaries 



The Yale code (iDemargue Sz Fercvl Il964l: 



Guenther. Demarque. Kim Sz Finsonneaultl 19921 ) 
uses as an outer boun dary condition the ein - 
pirically derived fit of Krishna Swamv ( 19661 ) 
to the T — T relation for the Sun, e Eridani, 
and Gmb 1830. Empirical fits have the flaw 
that they are suspect if extrapolated; these 



^Our goal is to develop a software environment that 
allows modification of physical modules by logical 
switches, thus maintaining consistency between old 
and new implementations. We plan to persist until 
we have an option that removes even the small differ- 
ences which remain for the standard solar model. 
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Table 3 

"Blob Sizes" and Mixing Length 



Model 


OiML 


9ml 


h/Hp 




A 


1.650 


1.0 


1.65 


1.0 


B 


2.323 


42.0 


0.358 


0.154 


C 


3.286 


270.0 


0.200 


0.0608 


D 


4.000 


595.0 


0.164 


0.0410 


E 


5.190 


1,540.0 


0.132 


0.0255 




log Te(K) 



Fig. 1. — Solar models which differ only by 
the mixing length and the geometric factor, 
scaled to have the same radius. The trajecto- 
ries in the HR diagram nearly overlay one an- 
other except on the Hayashi track, where they 
are poorly c onstrained. Note that the stel- 
lar birthline dStahler k Fallal l2004l ^ lies near 
logL/L0 ~ 1; the more luminous parts of the 
tracks ignore accretion and so are not realis- 
tic. 



stars are on the main sequence, and of G 
and K spectral type (G2V, K2V, and G8Vp, 
respectively). Gmb 1830 is a halo star of 
0.64M. P) with a metalicity of about 0.1 of solar 



(|Allende Frieto. et al. Ilioool ^. while e Eridani 



is a solar metalicity star of about O.SSMq. If 
applied to stars of the same stage of evolu- 
tion and the same abundance, such empirical 
boundary condidtions are at their best. Un- 
fortunately the "calibration" approach may 
hide mistakes in the assumed physics. 

3.2. Atmospheric Outer Boundaries 



The G arching code (ISchlattll Il996l') was 



modified (jSchlattl. Weiss Ludwiglll997l ') to 
use synthetic atmospheres fitted to the inte- 
rior solution at optical depth (r = 20). In 
addition a spatially varying mixing length 
was employed to reproduce the pressure- 
temperature stratifica tion calculated by 2D- 
hydr odynamic models (jFreytag. Ludwig. &: Steffen 
19961 ). This involved the interpolation be- 
tween an atmospheric value (Balmer-line 
fits gave = 0.5) and an interior value 
(ctin t = 1.7 to get the correct solar radius) ; 
see ijFrevtag. Ludwig. Steffenlll996l ) for de- 
tails. This approach can be extended with a 
library of hydrodynamic model atmospheres 
(and unlike the Yale approach, is not in prin- 
ciple limited to G stars). However, we find 
that our own 2D simulations, because of the 
pinning of vortices, do not mix material as 
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efficiently as 3D. For a given driving, 2D 
gives higher velocities to maintain the same 
convective luminosity (Asplund. et al. 200d : 
Meakin fc Arnettll2007a| ]bl). Further. 2D sim- 
ulations have a different turbulent cascade 
and damping than 3D, which is related to 
this velocity difference. These issues need to 
be dealt with in making contact between ac- 
tual convective velocities and observed line 
widths. 

3.3. The Subphotospheric Region 

How does changing the oml — 9ml pair 
affect the structure of the sub-photospheric 
layers? Figure [2] shows models A through 
E in the log pressure — log temperature 
plane. Th is may be directly compar e d wit h 



Fig. 1 in ISchlattl. Weiss Ludwil (|l997l l. 



Model A is almost identical to their curve 
labeled "Eddington-approximation" , which 
used radiative diffusion and MLT with con- 
ventional parameters (essentially the same as 
model A, oml ~ 1-7 and qml = !)• In 
contrast, model D, which also used the Ed- 
dington approximation but used MLT with 
0(ML = 4.0 and qj^il = 595.0, closely resem- 
bles their curves labeled "2D-hydro-model" 
and "ID-model-atmosphere", and models C 
and E are similar. It appears that the signifi- 
cant point is not the choice of radiative diffu- 
sion versus radiative transfer, but rather the 



treatment of convection (|VandenBerg. et al 



20081 ) . The Yale group get hydrodynamics by 
empirical fitting to hydrodynamic observed 
atmospheres, the Garching group get hydro- 
dynamics by a fit to their 2D hydrodynamic 
atmospheres, and we get hydrodynamics by 
analytic theory based on 3D simulations of 
convection. 




log P 



Fig. 2. — Structure of sub-photospheric re- 
gions with different choice of the aML — 
9ml pairs. Model A is the lowest curve, 
and models B, C, D and E are suc- 
cessively higher. Model A is similar 
to the "Eddington-approximat ion" case of 
Schlattl. Weiss fc Ludwi3 (|l997l ) while mod- 
els C, D, and E are similar to their "2D-hydro- 
model" case. 




Fig. 3. — Structure of the super-adiabatic re- 
gion (SAR), with AV = V — Vad versus log- 
arithm of pressure (dynes/cm^). The width 
of the SAR is much less than a pressure scale 
height; this may be compared with i^/Hp in 
Table 2. The small blob size implied in mod- 
els C, D, and E are consistent with the small 
width of the SAR, which is all we should ex- 
pect without a 3D atmospheric model. Notice 
the small sub-adiabatic region just below the 
photosphere (the left of the graph, indicated 
by the arrow). 
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4. The SAR and Surface Velocities 
4.1. The Geometric Factor qml 

Although the traditional procedure for cal- 
ibrating stellar convection is the variation of 
the parameter auL^o adjust the stellar radius 
keeping qml constant, this is not the most 
natural choice. It is gML that determines 
the radiative diffusion rate from "convective 
blobs" , and is most effective in the super adi- 
abatic region (SAR). In the adiabatic regions, 
MLT gives an adiabatic gradient, so the choice 
of aML is irrelevant to structure there. His- 
torically, the reasonable choice — of forcing 
a one-parameter family by assuming constant 
values for all parameters except uml — has 
obscurred the physics. Simulations uncovered 
this mistake, with the indication that oml is 
determined by the dissipation w hich is fed by 



the turbulen t cascade, exactly as iKolmogorov 
(Il94ll . Il962l ) suggested. 



The geometric factor gML may be ex- 
pressed in terms of a ratio of time to transit 
a mixing length to time for diffusion to re- 
move the super-adiabatic exces s from a " blob" 
( Kippenhahn &: Weigert 199d ). It is not well 
defined because of geometric vagueness about 
the "blob"; here we take it to be propor- 
tional to the inverse square of the ratio Ib/^-m-, 
where liy is the blob diameter and im the mix- 
ing length. This is a deviation from MLT, 
for which = i^- With this identification 
we can compare the blob sizes for different 
otML-QML combinations given in Table [2j 

This is shown in Table [3j Notice that 
for larger values of mixing length parameter 
aML-, the blob size becomes smaller, whether 
measured relative to a pressure scale height 
ih/Hp = 1/ ^qml or relative to a mixing 
length Ih/^m = ^/{aML^JguL)- This means 
that, for acceptable solar pairs of aML^OML, 
larger values of the mixing length imply nar- 
rower and more intense superadiabatic regions 
to drive the convection. Larger values of mix- 



ing length parameter aML imply larger veloc- 
ity scale (larger Vm) as Table [2] shows. Thus, 
models A-E are a sequence having increas- 
ingly vigorous and narrowly restricted regions 
of convective driving (acceleration). 

Figure [3] shows the structure of the SAR for 
models A t hrough E. This may be co i npare d 
to Fig. 2. of lSchlattl. Weiss fc Ludwi3 (|l997l l. 
Again model A resembles their "Eddington- 
approximation" curve, and models C, D and 
E are similar to their "2D-hydro-model" and 
"ID-model-atmosphere" curves. Here AV = 
V—Va is plotted against logarithm of pressure 
(dynes/cm^). 

Above the horizontal line AV = 0, buoy- 
ant forces accelerate the turbulent flow, while 
below the line we have buoyancy damping 
(deceleration; this region is barely visible at 
the left edge of the curve). According to 
MLT with the Schwarzschild criterion for con- 
vection, there should be no flow for AV < 
0. The area above (under) the curve gives 
the net buoyant acceleration (deceleration). 
Clearly the deceleration, seen as the small 
depression near logP = 5, is overcome by 
the much larger region of acceleration around 
logP ~ 5.2, so that the Schwarzschild crite- 
rion gives incorrect results here. The area ar- 
gument implied in Figure [3] is esser itially the 
bulk Richardson number crite rion ( Fernandol 
199ll : iMeakin fc Arnettll2007bl '). and is nonlo- 
cal. Therefore, because the pathological de- 
celeration implied by use of the Schwarzschild 
criterion is incorrect, the velocities Vm given 
in Table [2] are directly related to those which 
produce solar line broadening. 

4.2. Micro- and Macro-turbulence 

FigurelHshows the run of turbulent velocity 
as a function of optical depth for the five mod- 
els, and for the Nordlund-St ein 3D hydrody- 



namic atmosphere quoted in lAsplund. et al. 
Jiooi). The s emi-empirical st e llar a tmo- 



sphere models of Fontenla et al. ( 20061 ) give 
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Fig. 4. — Convective velocities versus log op- 
tical depth for solar models which differ only 
by the mixing length and the geometric factor. 
The convective velocity changes while there 
are no other significant changes; standard 
mixing length theory with the Schwarzschild 
criterion was used. Case A had ajv/L = 1-643 
and the usual geometric factor, gML = 1- 
Case D {oml = 4 and gML = 595.0) is the 
estimated value for saturation of the dissi- 
pation length for a solar convection zone of 
depth of 20 pressure scale heights. In MLT 
the velocity scale is not constrained physically, 
but only fixed by historical parameters (which 
are inconsistent with both 3D simulations and 
hydrodynamic th eory). The 3D niqdel a tmo- 

are 



sphere data from Asplund. et al. ( 20051 ) 



dramatically different at small optical depth. 



curves similar to those of Nordlund-Stein, but 
are not plotted to reduce crowding. The most 
striking feature in this figure is the differ- 
ence between the low depth behavior of the 
models (an abrupt cliff at r ~ 1), and the 
3D-atmospheres (a gentler slope for lower r). 
This is due to the use in the ID models of the 
Schwarzschild criterion for convection, a local 
condition. A weakly-stable stratification can- 
not really hold back vigorous motion, as use 
of the local Schwarzschild criterion implies. 

The micro- and macro-turbulent velocities, 
(mi and Crna, are parameters which were intro- 
duced long ago El to account for the embarass- 
ment that, according to the Schwarzschild 
criterion, conventional solar atmospheres are 
not co nvective a t the surface. Note that if 
C = JCL + Cla, then 1.9 < C < 3.0 km/s for 



the Sun (jcoxlbood ^. This is indicated by the 
vertical bounded line in Figure HI The con- 
nection between this C and the actual turbu- 
lent velocity due to convection is not simple, 
involving line-formation, photon escape, and 
inhomogeneous stellar surface layers. 

Fortunately, mu lti-dimeiisional hydrody- 

namic atmospheres ( AsplundlboOO ; Nordlund Sz Stein 
2003) do provide a spectacular fit to line 



shapes, with no free parameters, so we iden- 
tify the convective velocities well below the 
photosphere (optical depth r > 3) in these 
simulations with those predicted by our 
hydro-dynamically consistent choice of mix- 
ing length parameter oml- This means we 
are essentially matching different 3D simula- 
tions in the region of adiabatic convection, 
where they should give identical answers, and 
minimizing the sensitivity of the match to 
the complexities of atmospheric detail. Op- 
tical depth is sensitive to temperature (the 
opacity is k oc here), so that the visible 
surface is a complex structure (see Fig. 24 



^See iHuang fc Str"^ l|l96d ) for an early review, in 
which is already a well established parameter. 
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in Stein &: Nordlund ( 19981 )). For example, 
a 10% fluctuation in temperature implies a 
change in 2.4 in the opacity. The optical 
depth of the photosphere occurs at different 
radii for different positions on the solar sur- 
face, so that fitting it with a single radius 
is difficult. At greater depths we expect the 
3D atmospheres and the ID models to agree, 
but near the surface it is not clear that the 
3D and ID definiinitions of optical depth are 
consistent. 

Pressure should be a better coordinate for 
matching 3D results to a ID model. Unlike 
the optical depth, the pressure is a weaker 
function of angular position on the solar sur- 
face. Hydrodynamic flow tends to smooth 
pressure variations, making the deflnition of a 
mean pressure-radius relation more meaning- 
ful. Figure [5] plots convective velocities ver- 
sus log pressure for models A through E. We 
can see that the lAsplund. et al. I (l2005l ) model 
smoothly joins onto model D. 

Let us construct a simple model of the 
motion in this region to see how hydrody- 
namic arguments might give modiflcations to 
the purely hydrostatic boundary conditions 
used in models A through E. We will assume 
that the velocity is dominated by flow at the 
largest scales of turbulence. These scales con- 
tain most of the kinetic energy, and are least 
non-laminar. Convective motions are driven 
by the sinking of matter which is cooling due 
to transparency near the surface. This gener- 
ates gravity waves in the near-surface region. 
We will approximate the large scale average 
of this motion by g-modes (jLandau &: Lifshitz 
( 19591 ). see § 12) whose amplitude falls off ex- 
ponentially with pressure scale height. This 
implied a scaling with position above an in- 
terface at radius ro, P{r) = P{rQ)v{r) /v{rQ). 
Despite its extreme simplicity and harsh as- 
sumptions, this simple picture gives a signif- 
icantly improved approximation to the be- 
havior of the velocity in the photospheric re- 



> 1- 





1 1 1 1 1 1 1 1 . 


3D atmosphere / 




Simple hydro LJ 





4 4.5 5 5.5 6 

log P 

Fig. 5. — Convective velocities versus log 
pressure, for solar models and 3D hydrody- 
namic solar atmospheres. The atmospheres 
extend to lower pressure than the solar mod- 
els (actually the atmospheres extend to higher 
pressure too, but these va l ues w ere not in Ta- 



ble 1 of lAsplund. et al. I (|2005l )). It is clear 



that the 3D atmospheres would join smoothy 
onto solar model D for oa/l = 4.0 and qml = 
595.0, as we would have predicted. The thin 
solid curve labeled "Simple hydro" represents 
a hydrodynamic extrapolation from the point 
of maximum convective velocity (see text). 
Replacing the MLT estimate (based on the 
Schwarzschild criterion and hydrostatic struc- 
ture) with a physically motivated extimate 
gives a strikingly better agreement with both 
the 3D atmospheres and the empirical solar 
data. 
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gions. The thin black hne labeled "Simple 
hydro" in Figure [5] represents such flow, fit- 
ted from the point of maximum convective 
velocity in model D. It captures the quali- 
tative behavior far better than the conven- 
tional hydrostatic assumptions (shown as the 
steep "clifF' near logP = 5.25), and promises 
to do better as the complex physics of the 
photosphere is more fa it hfully represented 



(IStein Nordlund 
2000l ). 



19981 : iNordlund k Stein 



This suggests that the photospheric veloc- 
ity may be estimated by C ~ 0.8 Vm, which 
predicts a connection between fitted line 
shapes and convective flow. Better physics 
for turbulent flow seems to be needed in 
ID stellar atmospheres, and some 3D fea- 
tures are difficult to represent in ID, such 
as inhomogeneity be tween upward and down- 
ward moving fiow s (IStein &: Nordluiid Il998l: 



Nordlund fc SteinI l200d : ISteffen. M 



2003). 



The 3D hydrodynamic atmospheres can pro- 
vide insight into the correct mapping of real- 
istic physics of a multi-modal region onto a 
ID stellar model, and tighter constraints on C, 
for a given Vm- For Models A through E, this 
condition favors Model D. 

Independent of a ny estimate of C, our sim- 
ulations and theorv (iMeakin &: Arnettll2007bl : 



Arnett. Meakin. &: Young] |2009| ) suggest from 
hydrodynamics alone that models C, D and E 
are most plausible, i.e., oml lies in the range 
of 3 to 5 because of enhancement of turbu- 
lent damping in deeply stratified convection 
regions (a ml "saturation"). This consistency 
is encouraging. 

5. Summary 



simulations and theorv (Meakin &; Arnett 


2007b: 


Arnett. Meakin. & Young 


2009) have 



been applied to sub-photospheric regions of 
solar models. Even within MLT, a dynami- 



cally consistent velocity field (i.e., a consistent 
choice of oml and an adjusted qml), gives a 
better agreement with 

1. empirical T(r) relations, and 

2. 3D hydrodynamic models of stellar at- 
mospheres. 

Using the correct condition for mixing (the 
bulk Richardson number) implies that the 
ID atmospheres should exhibit hydrodynamic 
fiow. Fu rther, simpl e hydrodynamic consid - 
erations ( Press 1981 : Press &: Rvbicki 198ll ) 
suggest g-mode waves will be generated and 
penetrate to the photosphere (these are gen- 
erated by turbulent forcing from convection). 
We show that there is a connection between 
the predicted turbulent velocity scale and the 
observed (macro and micro)-turbulent veloc- 
ities, which removes the embarrassment of 
non-convective surface regions predicted by 
ID stellar atmosphere theory. As a bonus, 
we find that the observed macro- and micro- 
turbulence for the Sun can be used to fix the 
choice of qml (model D). 

We may also have a resolution of an ap- 
parent contrad iction. Atmos p heric models of 



white dwarfs (IWingetl l2008l : Bergeron et al 
19951 ). which have shallow convection zones, 
use MLT parameters (ML2: a = 0.6, con- 
siderably smaller than used for the Sun), 
indicating less vigorou s convection. Low 



2005 



mass eclipsing binar ies (jStassun. et al. 
Morales, et al. 20081 ) are generally fit with 
a ~ 1 (again low convective efliciency), these 
models do not have shallow convection zones. 
In MLT there is no rationale for these differ- 
ences. Use of Eq. [T] will give models having 
thin convection zones which agree with MLT 
models using small omLj so we expect to re- 
produce the white dwarf results. For low mass 
stars, the surface temperatures will be lower 
than the solar value, so that the SAR should 
comprise more mass, i.e., we expect larger 
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9ml to be physically correct. Table [2] indi- 
cates that there is a trade off between oml 
and qml- to compensate for lower gML, a 
must be lower, for the same radius. For a 
deep convection zone, oml is fixed; then a 
stronger SAR (larger gML, and more inef- 
ficient convection) will give a larger radius. 
This is the sense of the discrepancy of the 
computed radii for low mass eclipsing binarie s 



(|Stassun. et al.1 I2OO8I : [Morales, et aD boosl l 



and we suggest that part of the discrepancy 
may be due to the convection algorithm used. 
Unfortunately, direct calculation of low mass 
dwarfs (M ^ O.2M0) with oml ~ 4 exposes 
limitations in MLT: the SAR is forced upward 
into the photosphere, making 3D atmospheres 
a necessity for gaining insight into a plausible 
treatment in stellar models. 

We have, in fact, sketched a way to elimi- 
nate astronomical calibration from stellar con- 
vection theory: 

1. Adjust uml from convection simula- 
tions. The mixing length is im = 
OiML Hp (where Hp is the pressure 
scale height), and equal to the depth of 
the convection up to 4 Hp, and oml ~ 4 
for deeper convection zones. 

2. Adjust gML from 3D hydrodynamic at- 
mosphere simulations, fitting the curve 
of super-adiabatic excess in the supera- 
diabatic region. This is more accurate, 
but equivalent to adjusting qml to re- 
produce a self-consistent SAR. 

Notice that a fit to the present day solar ra- 
dius is not logically necessary. 

By seriously considering MLT, we have de- 
termined that no significant free parameters 
are left to adjust within the framework of the 
theory. We find that the choice of two charac- 
teristic lengths, which are determined by the 
flow, close the system: the turbulent dissipa- 
tion length and the size of the super-adiabatic 



region (SAR). Alternatively, the constraint 
that the observed micro- and macro-turbulent 
velocities agree with those predicted using the 
bulk Richardson criterion for surface convec- 
tive mixing can be used instead of the SAR 
size. 

MLT is still an incomplete theory, but it 
is suggestive that even modest changes to- 
ward a better physical interpretation, based 
upon 3D simulations and on a more complete 
turbulence theory, do give improvements in 
the models. MLT, as used here, may be de- 
rived from a more general turbulent kinetic 
energy equation by ignor i ng ce rtain terms 
(|Arnett. Meakin. &: Youn3 bOQ^ ). Some of 
the ignored terms are important, emphasiz- 
ing that MLT is incomplete. However, the 
approach sketched above may be generalized, 
and inclusion of missing terms gives a con- 
vection theory that is nonlocal, time depen- 
dent, provides robust velocity estimates, and 
is based on simulations and terrestrial experi- 
ment, with no astronomical calibration. This 
more difficult theory will be presented in de- 
tail in future publications. 

This work was supported in part by NSF 
Grant 0708871 and NASA Grant NNX08AH19G 
at the University of Arizona. We wish to 
thank Robert Stein for discussion of his un- 
published work on turbulent damping in solar 
convection simulations, Robert Buchler for 
discussions on modeling time-dependent con- 
vection, and Martin Asplund for providing 
machine-readable copies of solar surface mod- 
els. 
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A. MLT Geometric Parameter 

This analysis uses the formulation of Kippenhahn &: Weigert ( 19901 ): see their discussion for more 
detail. In the mixing-length theory, there are two important conditions which involve radiative 
diffusion: luminosity conservation and blob cooling. The simple condition L = L{rad) + L{conv) 
is written as 



(V- V, 



U{Vr-V), 



(Al) 



which is identical to Eq. 7.15 of Kippenhahn Sz Weigert ( 1990l ). Here the subscripts on the V's 
denote e for mass element (the blob), a for adiabatic, r for radiative, and no subscript for the 
background (environment) value. The diffusive cooling of the blob implies 



Va=gML2UiV -VeP, (A2) 

which is identical to Eq. 7.14 of Kippenhahn &: Weiger 

3 (|l99d ). except for the introduction of a 
scaling factor qml- For qml = 1 we regain conventional MLT. Thus, the definition of U becomes 



9mlU 



9ml 



3acT^ .8Hp\ 



(A3) 



which is their Eq. 7.12 with an extra factor qml-, and our j3T is their 5. If we define U* = qmlU 
and = V - Va + (C/*)^, we may write 



*\3 



+ 



^QML 



u*{c - {u*Y - VF) = 0, 



(A4) 



which is Eq. 7.18 of Kippenhahn &: Weigert ( 19901 ). except for the factor of Qml in the denominator 
and the replacement of U hy U* . The same solution procedures may now be applied to solve for C, 
and hence V. Any value of Qml that is not excessively large or small (within a few powers of ten 
of unity) has no significant effect except in regions that are both convective and nonadiabatic. 

An estimate of qml in terms of the size of a convective "element" or "blob" is given in Tabled] 
above, which we repeat here: qml = (^/v^^^fe)^, where is the "blob" radius. In MLT, = 
£/\/3 K, 0.577^, forcing two independent length scales, I and y/?>ri„ to be the same. 

Adjustment of qml allows the super-adiabatic region to have the correct entropy jump, for any 
reasonable value of the mixing length parameter aML'-, that is, uml may be chosen to be hydro- 
dynamically consistent. This does not provide a consistent convective theory if other important 
effects, such as entrainment and wave generation, are ignored. 
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